AlAA-98-2022 

FRACTURE MECHANICS ANALYSIS OF 
STITCHED STIFFENER-SKIN DEBONDING 

E.H. Glaessgen,* I.S. Raju f and C.C. Poe, Jr/ 

NASA Langley Research Center, Hampton, VA 23681-0001, U.S.A. 


Abstract 

An analysis based on plate finite elements and the 
virtual crack closure technique has been implemented to 
study the effect of stitching on mod^ I and mode II strain 
energy release rates for debond configurations. The 
stitches were modeled as discrete nonlinear fastener 
elements with a compliance determined by experiment. 
The axial and shear behavior of the stitches was 
considered, however, the two compliances and failure 
loads were assumed to be independent. Both a double 
cantilever beam (mode I) and a mixed mode skin- 
stiffener debond configuration were studied. In the 
double cantilever beam configurations, Gi began to 
decrease once the debond had grown beyond the first 
row of stitches and was reduced to zero for long 
debonds. In the mixed-mode skin-stiffener 
configurations, Gi showed a similar behavior as in the 
double cantilever beam configurations, however, Gn 
continued to increase with increasing debond length. 

Introduction 

Stitched warp-knit textile composite materials are 
currently being considered for use in primary aircraft 
structures. 1 In the NASA Advanced Subsonic 
Technology (AST) program, a stitched composite wing 
skin is being developed to demonstrate both the 
manufacturing and analytical technology needed to 
produce such structures. Structural panels, such as the 
one shown in Figure 1, are made of stiffening elements 
that are mechanically attached to the skin with stitches 
in the dry preform state before the resin is introduced via 
resin film infusion (RFI). 1 

The panel configuration shown in Figure 1 contains 
several important features: the skin, stiffeners, and 
intercostals. These panels are typically made of between 
two and ten stacks of 0.055 in. thick carbon warp-knit 
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fabric that are layered and stitched with Kevlar yams to 
form a skin. The stiffeners and intercostals are made of a 
similar number of stacks of stitched fabric and are 
stitched to the skin. Once the preform is assembled, the 
entire structure is infused with epoxy resin. 

Stiffened panels are typically subjected to large in-plane 
and out-of-plane loads that produce considerable bending 
and shearing stresses at skin-stiffener interfaces. 2 Much 
of the research on skin-stiffener debonding has focused 
on the calculation of these skin-stiffener interface 
stresses. The stresses may be large enough to cause a 
separation between the skin and stiffening elements 
resulting in a delamination or debond. Fracture 
mechanics-based approaches have been widely used to 
characterize delaminations and debonds in unstitched 
structures. Two-dimensional plane strain models have 
been used for analysis of skin-stiffener debonding. 6 ' 8 
Models based on quasi-3D or 3D brick finite element 
models have been used to study edge delamination and 
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near-surface delamination of composites. Since many 

layers of brick elements are often required to model 
both the skin panel and the associated stiffeners, finite 
element models with large numbers of degrees-of- 
ffeedom may be required. 

In an effort to develop computationally efficient models, 
references 12- 15 proposed the use of plate elements to 
model skin-stiffener debond problems and calculate 
strain energy release rate using the virtual crack closure 
technique (VCCT). Conventional plate modeling 
inherently assumes that the reference surface of the plate 
coincides with the middle surface. Thus, the skin and 
stiffener are usually modeled by plate elements with 
nodes at their respective mid-planes. This conventional 
method is not convenient for modeling debonding 
because it entails complex constraints to tie the flange 
nodes to the corresponding skin nodes. Thus, the 
approach taken in references 12-16 and the present 
analysis, is to place the skin nodes and the stiffener 
nodes along the interface between the skin and the 
stiffener. The positioning of these nodes at the interface 
is performed by defining an offset distance from the mid- 
plane of both the skin and the stiffener. Plate-element 
models using this technique can be used to evaluate 
accurate values of mode I and mode II strain energy 
release rates. 

Suppression of the growth of debonds between the 
stiffeners and skin may be achieved by stitching the 
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stiffeners to the skin. The effects of stitching on 
delamination or debond growth in composites have 
been examined in simple configurations with some 
success. In these works, the stitches are modeled as 
truss or beam elements connecting nodes through the 
thickness of the material. References 1 7-2 1 modeled 
laminates as two-dimensional plane strain structural 
components, while references 22-23 modeled the 
laminates as three-dimensional solids. An advantage 
of the three-dimensional modeling is to allow the 
stitches to be modeled discretely rather than as 
structural components with an “effective” stiffness. 

The objective of this paper is to Quantify the effect of 
stitches on the mode I and mode II strain energy release 
rates of double cantilever beam and mixed mode skin- 
stiffener debond configurations. The plate element 
modeling technique is used to analyze the debond 
configurations and the VCCT is used to calculate strain 
energy release rates. Because the nodes of the plate 
elements are offset to the interface between the skin and 
stiffener, the plate element models do not allow for 
nodal connectivity beyond the element interface along 
the debonds. The stitches are not modeled as finite 
length spar or beam elements as in references 17-23, but 
rather as nonlinear fastener elements with axial and shear 
compliances determined by experiment. In this paper, 
the effect of parameters such as stitch stiffness, applied 
load and debond length on the strain energy release rates 
are studied. 

Analysis 

In this section, the analyses for the double cantilever 
beam and debond configurations are presented. Next, 
the method used for computing strain energy release 
rates using the VCCT is briefly discussed. Lastly, 
strategies for modeling the stitches are presented. 

Double Cantilever Beam and 
Mixed-Mode Debond Configurations 

Figure 1 shows the configuration and loading of a 
typical skin-stiffened stitched composite panel. Three- 
dimensional modeling and analysis of this complex 
configuration would require a large finite element model 
with many thousands of degrees of freedom. However, 
considerable insight into the behavior of the complicated 
configurations can be gained by examining simpler 
configurations such as those shown in Figure 2 while 
reducing modeling complexity. Figure 2 shows the two 
simple configurations, the double cantilever beam 
(DCB) and the mixed-mode skin-stiffener debond 
configuration. Both of the configurations were modeled 
as infinitely wide strips with cylindrical bending 
repeating unit boundary conditions (v=0, 0*=O on 
y=±S y /2 in Figure 2). 


For both the DCB and mixed-mode debond 
configurations, a stitch pattern of 8 rows/in. by 8 
stitches/in. was assumed. That is, the stitch spacings, 
s x and Sy, are 0.125 inches in the x- and y- directions, 
respectively (see Figure 2). Axial and shear stitch 
stiffnesses and forces are denoted by K ax iai. K she ar and 
F axial andFshear, respectively, in Figure 2. Material and 
skin thicknesses that are representative of the stitched 
composite wing skin used in the NASA Advanced 
Subsonic Technology (AST) program are considered. 1 

A representation of the finite element mesh used in these 
analyses is shown for both the DCB and mixed-mode 
debond configurations in Figures 2(a) and 2(b), 
respectively. The element size chosen was 0.005 in. in 
the x-direction by 0.03125 in. in the y-direction for both 
models. 

The material is a carbon/epoxy system representative of 
materials used in the AST stitched composite wing 
skin. 1 In both the DCB and mixed-mode debond 
analyses, the same material is assumed for both the skin 
and the stiffener flange. Each stack of material is 
oriented with its primary axis in the x-direction and is 
composed of forty-four percent zero degree fibers, forty- 
four percent forty-five degree fibers and twelve percent 
ninety degree fibers. The equivalent laminate stacking 
sequence is (45/-45/0/90/0/-45/45)„ s with areal weights 
of 0.577 oz./ft. 2 , 1.21 oz ./ft. 2 and 0.651 oz./ft. 2 for the 
forty-five, zero and ninety degree plies, respectively. In 
these analyses, the laminates are assumed to be 
homogeneous with properties 

E„=9.25xlOV' E 3 3=E 2 2=4.67xl0V' 

|0.i 2 = pi3=2.27xl0 6 /?s/ |i.2 3 = 0.497xl0 6 /w/ 

Ui2= '1)13=0.397 -u 2 3=0.490 

where Ey, Py, Oy (ij=l,2,3) are the Young’s moduli, 
shear moduli, and Poisson’s ratio, respectively, and the 
subscripts 1,2,3 represent the fiber and two transverse 
directions, respectively. 

Strain Energy Release Rates 

24-25 

The virtual crack closure technique (VCCT) can be 
used to calculate strain energy release rates, G, with 
plate elements using the techniques discussed in 
references 14-16. The STAGS 480, 9-node quadratic 
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shear deformable, plate/shell element is used for 
modeling the debond configurations. A representation 
of the 9-noded plate elements near a debond front is 
shown in Figure 3. Reference 16 suggests that allowing 
the elements ahead of the debond front to have free 
rotations ensures accurate evaluation of the strain energy 
release rates. If there are free rotations ahead of the 
debond front, then the G values can be calculated using 
the nodal forces (F x , F y , F z ) and displacements (u, v, w) 
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near the debond front and the increment of new debond 

12 

area created as (see Figure 3) 

Mode-! components: 


(G[)i = -■ 
«?/); = - 


2Ab~[ Fz ‘( Wp ~ W P') + F zr(wi — W/')] 
^^-[ F zj{w q -w q )+ F lf (w m Wrri )]<» 

rrrtfjw, - wv) + F z ,(w„ ~ w„)] 
It-sbk 


Mode-U components: 


this modeling strategy, that assumes no rotational 
constraints ahead of the debond front, is termed 
“Technique-B” in references 14-15. 

Both the models of the double cantilever beam and the 
mixed-mode skin-stiffener debond configuration assume 
self-similar debond growth between the skin and 
stiffener. Thus, no variation in G across the width of 
the model is assumed and the values of G reported are 
those calculated along the lines of the stitching (y = 0 in 
Figure 2). Although a tendency for the debond to leave 
the skin-stiffener interface has been shown in some 
composite skin-stiffener structures, 6 ' 8 experimental 
evidence suggests that the self-similar assumptions are 
valid for this type of stitched woven composite. 28 


(G„)i = --—[f x {u p - u P ) + F x ,{ui~ ui)\ 
l&bi 

(G„)j = -~^—[f XJ { u q - u <{ ) + F Xf {um - U m )\ (2) 
lusbj 

( G„)k = - [F xt ( Ur - Ur ) + F x ,(un ~ Un‘ )] 

ICsbk 

Mode-Ill components: 

(G m )i = -J~[Fy,(v P - v P ) + Fy t (yi ~ vr)] 

(G w )j = ~ v ?') + *J/(v- “ v m )\ ( 3 * ) 

(Gm)k = ~^^’[ F yS V r ~ Vr ) + Fy s {vn ~ V„)] 


with 

(Grotai)y = (Gi + Gii + Gut ) I r , and (4) 

Y=i, j and k. 


The elements are assumed to have the same length, A, 
ahead of and behind the debond front. The equivalent 
widths apportioned to the two comer debond-front nodes 
are b, and bk, and to the midside debond front node is 
bj. These are given by equation (5) as 


bi = -[bj-\ + bj], 
0 



bk = ^[bj+bj+\], 

0 


where bj. h bj and bj+i are the widths of layers J-l, J and 
J+l, respectively, as shown in Figure 3(b). Note that 


The strain energy release rates along the debond front of 
the DCB and mixed-mode skin-stiffener debond 
configurations are calculated using equations (l)-(5). 
The nonlinear nature of the compliance curves 
necessitates a geometrically nonlinear finite element 
analysis. Both of the configurations have been analyzed 
using the STAGS finite element code. 

Stitch Modeling 

Unlike the two-dimensional and three-dimensional 
models considered in references 17-21, the plate 
element-based modeling technique does not allow 
through-the-thickness modeling of details such as the 
stitches; nor does it allow nodal connections other than 
at the plate element reference surface. Thus, in the 
present technique, the stitches are not modeled as spar 
or beam elements, but rather as “fastener elements.” 
The fastener elements are imposed as nonlinear 
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constraints within the plate element model. These 
“fastener elements” have both an axial and a shear 
stiffness, Kaxiai and K S h C ar, respectively, and are 
schematically shown as springs in Figure 2(a). Only 
the fastener elements behind the debond front 
( L<x<L+a in Figure 2a) carry load since the upper and 
lower plate elements ahead of the debond front are 
coupled using constraint equations to have identical 
translational displacements. 

Accurate compliance curves for both axial and shear 
behavior of the stitches were developed in reference29 
using flatwise tension and double lap shear tests (see 
Figure 4), respectively. A piecewise linear 
representation of this data is used in the finite element 
model. The points used in the linearization of the 
compliance curves are shown in Figure 4. Examination 
of the flatwise tensile test data revealed that the axial 
compliance of the stitches increases from approximately 
1.30E-5 in./lb. (77000 lb./in. stiffness) initially to 
4.88E-4 in./lb. (2050 lb./in. stiffness) near failure. 
Similar examination of the double lap shear test data 
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showed that the shear compliance of the stitches 
increases from approximately 5.50E-5 in./lb. (18200 
Ib./in. stiffness) initially to 6.25E-4 in./lb. (1600 lb./in. 
stiffness) near failure. The axial and shear responses of 
the stitches are assumed to be independent. Also note 
that failure of the stitches occurred at 58 lb. in tension 
and 38 lb. in shear. These stiffnesses and failure loads 
will be used for the characterizations in this paper. For 
the purposes of these analyses, the axial and shear 
responses of the stitches are assumed to be independent. 

Results and Discussion 
Double Cantilever Beam (DCB) Models 

Effect of Stitch Stiffness 

A mode I loading of 125 lb./in. .was applied to the 
configuration shown in Figure 5(a) with a debond 
length, a, of 0.250 in. The axial stiffness of the stitches 
was increased in decades from 1.0 lb./in. to 100,000 
lb./in. to evaluate its effect on the DCB configuration 
with several sublaminate thicknesses, t, of 0.1 10, 0.220 
and 0.440 inches. Figure 5 presents the strain energy 
release rate, G, and the axial stitch force, FaxiJFpum, for 
the DCB configuration. A stitch stiffness in excess of 

1 .000 lb./in. is required to noticeably effect the value of 
Ci for even the most compliant (f=0.1 10 in.) DCB. As 
shown in Figure 5(a), for a given applied load, P, the 
effect of the stitching on Gi appears to decrease with 
increasing DCB thickness for the range of stitch 
stiffnesses considered. At stitch stiffnesses approaching 

100.000 lb./in. the strain energy release rate of the 
thinnest (most compliant) DCB approaches that of the 
stiffestDCB since the stitches have more effect in the 
compliant DCB configurations than in the stiffer 
configurations. As a result, the load carried by the 
stitches increases with increasing DCB compliance as 
shown in Figure 5(b). Comparison of the measured 
initial axial stitch stiffness of about 80,000 lb./in., from 
the flatwise tensile test data, with the curves of Figure 
5(a) reveals that this initial stitch stiffness will have a 
significant effect on strain energy release rate. 


Effect of Increasing Load 

Finite element models utilizing the compliance curves 
shown in Figure 4(a) were used to evaluate the response 
of the DCB to increasing load. As shown in Figure 
6(a), Gi increases dramatically with increasing applied 
load for all three DCB configurations. As expected, the 
thinnest (f=0. 1 10 in.) DCB shows the largest 
sensitivity to the increase in load. Figure 6(b) shows 


the increase in load, FaxJFjMum, in the stitches with 
increasing applied load for the three DCB 
configurations. However, Figure 4(a) shows that as the 
load in the stitch increases, the corresponding stitch 
stiffness decreases. Also note that the maximum values 
of strain energy release rate shown in Figure 6(a) are 
well beyond the value required to propagate the debond 
in an epoxy resin; Gic for a brittle epoxy resin is 
approximately 1.0 in.-lb./in. 230 

The effect of the nonlinearity of the stitches is most 
evident when the load in the stitches is allowed to 
increase to the large values shown in Figure 6(b). For 
example, in the thinnest (r=0. 110 in.) case with an 
applied load of 1,000 lb./in., the load of approximately 
42 lb. (FaxJFfaitur^ 0.73 in Figure 6b) in the stitches 
corresponds to a stitch stiffness of 2,600 lb./in. This 
represents a decrease of stitch stiffness by a factor of 30 
compared with the initial value. Thus, the effectiveness 
of each stitch is expected to decrease with increasing 
load. In other words, while the stitches significantly 
retard debond growth when they are carrying low loads, 
they have a decreased effecton strain energy release rate 
once they, and the debond front, are heavily loaded. 

Effect of Increasing Debond Length 

The third variable studied was the debond length. In 
this section, the debond length was allowed to increase 
with a fixed load of 500 lb./in. applied to the DCB. 
For the purposes of this discussion, a configuration 
corresponding to an intermediate DCB stiffness (/= 0.220 
in.) was chosen. As shown in Figure 7(a), an 
unstitched DCB exhibits an increasing Gi with 
increasing debond length, a. In contrast, the slope of 
the Gi curve for the stitched DCB begins to decrease 
immediately after the debond front passes the location of 
the first stitch (first vertical line in Figure 7a) and the Gi 
curve may actually begin to decrease for debond lengths 
longer than approximately 0.25 in. (i.e. a/t= 1.15). The 
decrease in slope of the Gi curve is larger with both the 
number of loaded stitches and the distance from the 
debond front to each of the load-bearing stitches. 
Further, for long debonds ( a/t>2.5 ), the strain energy 
release rate is nearly zero. 

The decrease in Gi for long debond lengths is the result 
of the stitches resisting debond opening (reducing the 
opening force in the region near the debond front.) 
This effect increases with increasing distance between 
the debond front and the stitch, and also with an 
increase in the number of load bearing stitches. Figure 
7(b) shows how each of the stitches begins to carry load 
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as the debond of length, a , passes its location in the 
model. For example, the force in the first stitch (Fmai, 
Stitch 1 in Figure 7b) starts at 0.0625 in. (a/t= 0.284) 
(point B) beyond the Stitch 1 location (point A). Since 
the results were evaluated from the finite element model 
with increments of debond growth of 0.0625in. 
(a/r= 0.284), the force in the stitch for debond lengths 
between points A and B is not known and is represented 
by the dashed line. The forces corresponding to the first 
0.0625 in. (a//=0.284) of debond growth beyond the 
location of each successive stitch are represented in a 
similar manner. As seen in the figure, for short debond 
lengths, only the first row of stitches is loaded. As the 
debond length increases, the load is shared by additional 
rows of stitches. 

Mixed Mode Debond Models 

The skin-stiffener configuration shown in Figure 2(b) 
exhibits both mode I and mode II deformations at the 
debond front. The configuration has a skin and a 
stiffener flange of equal thickness, t, of 0.220 inches. 
This mixed mode configuration is included to study the 
effect of increasing debond length on mixed-mode strain 
energy release rate and load carried in the stitches for a 
fixed applied load of 200 lb./in. As shown in Figure 
8(a), both the mode I and mode II strain energy release 
rates, G ( and Gn, respectively, for the unstitched 
configurations increase with increasing debond length 
over the range of debond lengths, a, considered. The 
effect of the stitching on Gi in the mixed-mode 
configuration is very similar to that shown for the DCB 
in Figure 7(a). Gi initially increases with debond length 
and then begins to decrease after the debond passes the 
first stitch (first vertical line in Figure 8a), decreasing as 
additional stitches begin to carry load, eventually 
reaching near-zero values for long debonds (a/t>2). In 
contrast, Gn increases slightly with increasing debond 
length over the range of debond lengths considered. 
Thus, even though G\ is significantly reduced, there 
may be sufficient Gn present to grow the debond. 

As seen in Figure 8(b), the axial force in the stitches 
( Fax, ad >n the mixed mode configuration shows a similar 
behavior with increasing debond length as in the DCB 
configuration. In this configuration, however, a shearing 
force (Fshear) is also present in the stitches, though it is 
smaller in both absolute and normalized magnitude than 
the axial force. As in the previous section, results were 
evaluated from the finite element model with increments 
of debond growth of 0.0625 in. (a/t= 0.284), so the force 
corresponding to the first 0.0625 in. (a/t= 0.284) of 
debond growth beyond a given stitch location in the 
finite element model was not recovered in the analysis 


and is represented by the dashed lines in the figure. 
There are two ways that the debond may continue to 
grow in this mixed-mode configuration. If the stitches 
remain intact, the debond may continue to grow by 
mode II as shown in Figure 8(a). However, if the 
stitches fail, non-zero mode I may be present and may 
also contribute to the growth of the debond. 

Concluding Remarks 

The effect of stitches on the mode I and mode II strain 
energy release rates of double cantilever beam and 
mixed-mode skin-stiffener debond configurations with an 
initial debond was studied. The plate element 
modeling technique was used to model the configuration 
and the virtual crack closure technique (VCCT) was 
used to calculate the strain energy release rates. Using 
this fracture mechanics approach, a debond between the 
skin and stiffener flange was assumed. The debond 
growth between the flange and the skin was assumed to 
be self-similar and continuous along the length of the 
flange-skin interface. The stitches were modeled as 
discrete nonlinear fastener elements with their 
compliance determined by experiment. Both axial and 
shear behavior of the stitches was considered, however, 
the two compliances and failure loads were assumed to 
be independent. 

Stitches with stiffnesses in excess of 1 ,000 lb./in. in the 
8x8 stitching pattern were required to noticeably affect 
the values of mode I strain energy release rate, Gi, for the 
double cantilever beam configurations studied. The 
effect of stitching was most significant for compliant 
DCB configurations. For fixed debond lengths, an 
increase in applied load results in an increase in the load 
carried by the stitches, thus, stitch compliance increases 
and contributes to an increase in Gi. However, if the 
debond length increases, additional stitches begin to 
carry load resulting in a decrease in G\. For long 
debonds, the stitches may produce enough compressive 
force at the debond front to close the debond and reduce 
Gi to zero. Thus, failure of the stitches may be required 
to propagate the debond. 

The mode I strain energy release rate computed for the 
mixed-mode debond configuration shows similar 
dependence on debond length as in the DCB. The 
stitches begin to carry considerable load and cause Gi to 
decrease once the debond length becomes sufficiently 
long. In contrast, the stitches have less effect on mode 
II and Gu increases throughout the range of debond 
lengths considered. Thus, there are two ways that the 
debond may continue to grow in this mixed-mode 
configuration. If the stitches do not fail, the debond 
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may continue to grow by mode II. However, if the 
stitches fail, non-zero mode I may be present and may 
also contribute to the growth of the debond. 
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Figure 1 Stitched composite panel with stiffeners 
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Dimensions not to scale 



Dimensions not to scale K ax j a |, K S h ear 



(b) Mixed-mode configuration 


Figure 2 Debond configurations 
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Tension force per stitch, lb. 

(a) Compliance of stitches in tension 


Failure Load=38 lb. 
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film 


8 stacks 


>(15, 0.0037) 


'(10, 0.00055) 15, 0.0017) 


' 3 stacks 


Shear force per stitch, lb. 

(b) Compliance of stitches in shear 

Figure 4 Stitch compliance (experimental results taken from reference 25) 
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Stitch Stiffness, K axja | (Jb./in.) 

(a) Strain energy release rate for a range of stitch stiffnesses 



Stitch Stiffness, K axia) (Ib./in.) 

(b) Stitch force for a range of stitch stiffnesses 

Figure 5 Effect of stitch stiffness on strain energy release rate and stitch force 
(Applied Load, p=125.0 lb./in., a=0.250 in., s x =0.125 in., s y =0.125 in.) 
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(a) Strain energy release rate as a function of applied load 
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(b) Stitch force as a function of applied load 

Figure 6 Effect of applied load on strain energy release rate and stitch force 
(a=0.250 in., s x =0.125 in., s y =0.125 in.) 
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Debond Length, a/t 

(a) Strain energy release rate as a function of debond length 



Debond Length, a/t 

(b) Stitch force as a function of debond length 

Figure 7 Effect of debond length on strain energy release rate and stitch force 
(Applied Load, p=500.0 IbVin., t=0.220 in., s x =0.125 In., s y =0.125 in.) 
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Debond Length, a/t 


(a) Strain energy release rate as a function of debond length 



(b) Stitch force as a function of debond length 
Figure 8 Effect of multiple stitches on G j, G|| and stitch force 

(Applied Load, p=200.0 IbTin., t=0.220 in., s x =0.125 in., s y =0.125 in.) 
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